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Abstract 



We consider here the structure of rotating compact objects endowed with 
a magnetic field in general relativity as models of pulsars. We discuss first the 
structure of rotating stars in the framework of Hartle taking different realistic 
equations of state and study their effects on bulk properties of the star. We 
consider the possibilty of rotating stars with a quark matter core. We further 
analyse the structure of the magnetic field in the interior of the star as affected 
by different equations of state as well as due to rotation. 
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I. INTRODUCTION 



The study of structure and dynamics of neutron stars is an interesting and challenging 
theoretical problem. The system here is quite complicated. The large density in the interior 
of such stars can lead to a possibility of melting the baryons to form a deconfined quark 
core IjT]]. Observations of rotating neutron stars with a period in the millisecond range has 
renewed the interest in the study of the influence and implications of rotation on the bulk 
properties of such objects [Q]. Further the neutron stars are also known to possess strong 
magnetic fields. Their polar field strengths are deduced from the observed spin slowdown 
of pulsars via the magnetic dipole breaking model with most of the young pulsars having 
a surface field, B ^ 10^^ gauss @]. Each one of the above facets of neutron stars poses 
interesting theoretical challenges. The fascinating possibility of a deconfined phase of matter 
in the interior of neutron stars has stimulated the work of many authors proposing different 
signals of existence of such a core of quark matter. These studies mostly have been done 
for considering nonrotating stars. The treatment of massive compact rotating objects in 
the framework of general theory of relativity is a rather complicated task. However, for 
slow rotation this can be done in the framework of Hartle and Thorne [Q. Further, the 
magnetic field structure outside the neutron star has been studied in connection with pulsar 
emission mechnism In a nonrelativistic approach, the generation and evolution of internal 
magnetic fields have been studied recently 0. A relativistic treatment to include structure 
of magnetic field in the interior of neutron stars by solving coupled Einstein Maxwell's 
equations has been done numerically in ref. 0. It is therefore interesting to incorporate all 
these effects together in its totality to the extent possible and the present work is a step in 
that direction. 

We shall consider rotating stars endowed with magnetic fields. The effects of rotation is 
taken into account in the framework of Hartle and Thorne where rotation effect is treated 
as a perturbation over the Schwartzchild solution and is valid naturally for slowly rotating 
stars The field structure shall be assumed to be poloidal and axisymmetric. The 
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effect of rotation of the star on the magnetic field structure is included, but not the reverse, 
amounting to the fact that the magnetic field energy is not strong enough to affect the space 
time geometry. We shall also take into account the effects of different equations of state on 
the gross properties of the star and on the magnetic field in the interior. 

We organise the paper as follows. In section II, we give a short description of treating 
rotating compact objects in general relativity in the approximations of Hartle and Thorne [§. 
In section III, we discuss Maxwell's equations in curved space time with Hartle metric and 
set up the differential equations for electromagnetic fields. Here we obtain the boundary 
conditions on the dipole magnetic field structure taking corrections due to rotation upto 
0(f2^). Using these boundary conditions, we integrate them from the boundary to the 
interior of the star. In section IV, we discuss different equations of state and their effects 
on the bulk properties of the star along with the field structure in the interior. Section V 
summarizes the results and discussions. 

II. ROTATING NEUTRON STARS IN HARTLE THORNE APPROXIMATION 

The description of rotating compact objects in the framework of general relativity is a 
complicated task. In what follows, we shall be using the method of Hartle and Thorne P|, 
appropriate for slowly rotating objects. Here a perturbative solution is formulated over the 
Schwarzschild metric for the nonrotating stars given as 



The metric function is expanded in terms of the star's rotational frequency, fi, the angular 
velocity of the star's surface relative to an observer at infinity and is given as 



In the above uj is the angular velocity of the star's fluid in a local inertial frame and depends 
on the radial coordinate r. The difference a) = f2 — u; is the angular velocity with which the 



ds^ = -e''"dt^ + (1 - 2GM/r)~^dr^ + r^{de^ + sin^edcj)'^) 



(1) 



ds' = -e^-dt" + e^^dr^ + e^^dO^ + e^^^dcj) - udtf 



(2) 
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fluid inside the star moves relative to a distant observer. The metric functions are expanded 
in multipoles around the Schwarzchild metric as 

e2-^e^°[l + 2(/io + /i2P2)], (3a) 

e^^ = (1 - 2GM/T)-' \l + — (mo + m^P^)/ ( 1 - ^"jl , (3b) 

L r V r /J 

e^'^ = [1 + 2{V2 - h2)P2] , (3c) 

e^^ = sin^ 9 [1 + 2{v2 - /i2)P2] • (3d) 

In the above,the perturbation functions mo, m2, Hq, /i2 and V2 are to be calculated from 
Einstein's field equations. Under the assumption that due to rotation the star distorts and 
affects the energy and pressure of the rotating fluid by Ae and Ap respectively i.e. the source 
term of Einstein equation — SttGT^,/ changes by AT^,^, one can perform a multipole 
expansion of these quantities as 

Ap= (e+p)[po+P2^2] (4a) 

Ae = (e + p)^bo+P;^2] (4b) 

Finally, the quantities ^'o(r), M(r) denoting the metric function and the mass profile of 
the nonrotating star, are obtained by solving the Tolman Oppenheimmer Volkhoff (TOV) 
equations 

dp G[e + p] [M + A'KT^p] 



dr r3 (r - 2GM) 
dM 



(5a) 



= A-nr e (5b) 



dr 

dvQ 2 dp 

dr (e + p) dr 



(5c) 
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In what follows, we shall be limiting rotational deformations to monopole approximation 
only. The differential equations for the perturbation functions are given as 

/ — \ 2 



4vrr2^(e + p)p* + -^j {^] - -^r=^^u;^ (6a) 



where, 



dr dp' "^"^^ UG-" \dr ) SG dr 

dpi _ Gmo{l + SnGr'^p) ATcGp^e + p)r 

1 r^f fdujV Id/ r^fuj^ \ 

^ 12(1-^^ {d^j ^3d^[i^nm)^ ^^^^ 



r 

Further, the expression for the angular velocity for the fluid relative to the local inertial 
frame, u is obtained from the equation 

1 d f A.duj\ idj „ , , 

r + = (6c) 

dr \ dr J r dr 

The metric perturbation function Hq is calculated using the algebraic relation obtained 
from hydrodynamic equation given as 

ho = -Po + Ve-'^'cu^ + hoc. (7a) 

Outside the star, mo and ho have the forms 

C P 

mo = SM--^, (7b) 

_ GSM GV^ 

°~ r-2GM r\r-2GMy ^ ' 

where, J is the angular momentum of the star given as 

The constant hoc in Eq. (|7aD is obtained by matching the solution for ho at the boundary. 
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III. ELECTROMAGNETIC FIELD STRUCTURE IN THE INTERIOR OF THE 

STAR 



Compatible with an axisymmetric and stationary system, one can write down the 
Maxwell's equations in a curved background specified by Hartle metric as given in Eq.(^). 
The nonzero components of the electromagnetic fields satisfy the equations P] 



_9 

dr 



+ 



d_ r 

de 



^u+X+tp+fj, jt 



(8a) 



d_ 

dr 



ijjt 



+ 



d_ r 



d_ 

dr 



+ 



d_ r 

de 



_gi.+A+V-+MJ</'^ (8b) 



d_ 

dr 



d_ 

^ de 



0, 



(8c) 



d_ 

dr 



d_ 

^ de 



e'^^Er 



d_ 

dr 



d_ 

de 



(8d) 



To consider the interior electromagnetic field structure, we shall assume the interior of 
the star to be a perfect conductor. Hence the electric field inside the star is zero and, the 
magnetic field components as modified due to rotation upto monopole approximation for 
the metric are given as 
d 



dr 
d_ 

de 



Ko/2^1 + /lo)5e] 







(9a) 



(9b) 



Assuming the magnetic field to be dipolar in nature, the field components can be written as 



Br = f{r) cose Bg = g{r) sine, 



(10) 



and the equations (^) and reduce to the coupled ordinary differential equations given 



as 



r(l + ho)g'{r) + 
r ) 



+ ( 1 1 il + ho + 



1 + ho) (l + ^z.^(r))+r/i[,(r) 
) f{r) = 



r - 



(11a) 



r/'(r) + 2/ + 2(l-^)"'^'(l 



Gmo 



r-2GM 



g{r) = 



(lib) 



The boundary conditions for these two coupled differential equations are found by match- 
ing the interior solutions with the exterior solutions at the boundary and are given as P,p!0 



f{R) = UR) + 5f{R), 



(12a) 



g{R) = go{R) + 6g{R), 



(12b) 



where, /o(-R) and 5'o(-R) are the solutions corresponding to the nonrotating star [T^ given 



as 



h{R) 



3fi 



In 1- 



2GMo 
R 



2GMo ^ GMo 



R 



R 



(13a) 



9oiR) 



3fi 



AG'^M^R 



1-^) ^ (13b) 



R J ' GMo""V R 

with, Mq = M{R) is the mass of the nonrotating star obtained by solving the TOV equations 
(pa|) and (|^). We now substitute equations (0) and (^3|) in equation ( p^T]) and expand in 
powers of GMq/R retaining terms upto monopole perturbation. Comparing the coefficients, 
we obtain the corrections to the nonrotating solutions in the monopole approximation as 

12/iG 



SfiR) 



i?4 



„ , , 1 4 GMo 

6M {- + 

4 5 R 



GJ2 / 1 3 GMo 

+ 



Vl4 10 R 



(14a) 



(5M 2 + 



37 GMo 

y R 



GJ^ /8 163 GMo 



i?3 



7^ 35 R 



(14b) 
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IV. EQUATIONS OF STATE AND QUARK HADRON PHASE TRANSITION 



The interior of the star has density which is much higher than the normal nuclear den- 
sities. Unfortunately no single description of the equation of state will suffice to describe 
the whole density range from the center to the surface of the star. The main problem in 
determining the equations of state (EOS) at high densities lies in the uncertainity of the in- 
teractions operative and the procedure to calculate the same. The scope of the present work 
will be limited to take different equations of state and their effects on different stellar prop- 
erties. We use here five neutron matter EOS and the bag model EOS for quark matter. The 
neutron matter EOS are (i) Wiringa's nonrelativistic EOS including the three body forces 
||Tl| ; (ii) Walecka's mean field theory EOS obtained with nucleons interacting with scalar 
and vector bosons |12|; (iii) a field theoretical EOS given for high density matter given by 
Sahu et al, assuming the composition to be neutron rich matter in beta equilibrium based on 
the chiral sigma model |]13|; (iv) Bethe- Johnson EOS based on a nonrelativistic constrained 
variational principle with a modified Reid soft core potential |jT^; and (v) Pandharipande's 
hyperonic matter EOS based on many body variational theory [|15 . 



For the densities below the nuclear density, the stellar matter is no longer determined 
by the nucleon nucleon interaction only and we shall use for the subnuclear density region, 
the EOS given by Negele and Vautherin for densities of 10^^ to 5 x 10^'^ gm cm~^, of 
Baym, Pethick and Sutherland [|l7l for densities of 5 x 10^° gm cm~'^ till ~ 10^ gm cm~^ 



and of Feynman et al [jT8| for densities below 10^ gm cm~'^. In Fig. 1, we plot the equations 
of state for different models that we shall be using to study the neutron star structure. An 
EOS is said to be stiffer if the pressure is higher for a given energy density. However, as 
can be seen from fig. 1 that none of the EOS is stiffer as compared to another in the whole 
range of the energy density considered. Thus the stiffness of an EOS can be defined in an 
averaged sense from the value of the maximum stable mass of the neutron star using the 
particular EOS. An EOS is said to be stiffer on an average if it has higher value for Mmax- 
In our calculations, however, we shall be referring to the stiff or soft nature of the EOS 
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corresponding to the central energy density chosen for the star. 



A. Hybrid stars 

Here we shall consider the case of a hybrid star i.e. a compact star with a quark matter 
core and a neutron matter crust. We shall be following the simplest approach of Ref ||T9( 



where one uses Bethe- Johnson EOS for neutron matter and bag model EOS for the quark 
matter phase. Bethe- Johnson EOS is based on a nonrelativistic constrained variational 
principle with a modified 'Reid' soft core potential ||2^ as input and reads ||21|| : 



e(n) = n[940.6 + 247.5n^-^^] p{n) = 366.3n^-^\ (15) 

where e{n) is the energy density, p{n) the pressure in MeV/fm^ as a function of baryon 
number density, n in fm~'^. The quark gluon plasma phase is taken into account through a 
bag model EOS given as 

p =1(^-45), (16) 

with the energy density of the QGP as that of a mixture of gluons and massless u-, d- and 
s- quarks at zero temperature given as 

eM = ^{l--as)fi% + B. (17) 

In the above, fis is the baryon chemical potential, is the strong coupling constant and 
B is the bag constant. As in Ref. [1^, we shall also take the strange quark mass to be 



zero. A finite strange mass will slightly decrease the pressure, which is inappreciable as 
compared to the role played by B and as ||T9| , p2| . Further, since we take the mass of the 
strange quarks to be zero, the matter is electrically neutral without any further need for 
electrons ||23|. Since we have only one chemical potential fiB, the phase boundary between 
hadronic and the QGP phase is determined by the Gibbs criteria as given by 

Phadron ifJ'cr) = PggpifJ'cr) (18) 



where Hcr is the critical baryon chemical potential. 

In Fig. 2, we consider such a phase transition construction for the bag model EOS for 
quark matter and the Bethe- Johnson EOS for the neutron matter. We have taken the bag 
constant as 5^/^=190 MeV and the strong coupling constant 0^=0.4. This leads to a critical 
chemical potential to be 1600 MeV and critical pressure, 407 MeV/fm^. A reduction of B 
leads to a decrease in the critical chemical potential. 



V. RESULTS AND DISCUSSIONS 



In order to construct a stellar model we proceed as follows. With a given central energy 



density as input, we integrate the coupled TOV equations (|5aD and (^) radially outwards 
from the center, using an equation of state until the pressure vanishes which defines the 
surface of the nonrotating star. This gives the mass (Mo) and radius {R) of the star as well 
as the pressure and mass profiles (p(r), M(r) respectively) to be used for later computations. 
The metric function z/o(r) inside the star is then calculated by integrating Eq. ( pcf ) from the 
surface to the center with the initial condition as 

2GMn 



uo{R) = In (^1 



R 

We next proceed to calculate the rotational perturbations by integrating equations 



(pbf), ( pq) from center to the surface for a given angular velocity at the center Uc with the 
initial conditions 0] 

^|.=o = 0, mo(0) = 0=p*(0), 

using the mass and pressure profiles obtained from solving the TOV equations. The second 
initial condition obviously implies that the central energy density of the rotating star is same 
as that of the non-rotating star [^. 

With the monopole corrections calculated as above, we then proceed to integrate the 
dipole magnetic field configuration functions /(r) and g{r) given in eq.s (|llb| ) and (|lla]) 
along with the respective boundary conditions given in Eq. ( |12a| ) and ( |1 2b|) . 
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The results with different equations of states are now given as follows. 

We have plotted in fig. 3a the energy density profiles for stars with different neutron mat- 
ter equations of state choosing the central energy density of the star as ec=1600 MeV/fm^. 
The drop in the profile near the surface corresponds to the crust of the star. The thickness 
of the crust increases with the softness of the EOS used to study the star The energy 



density profile is plotted for a pure quark star and a star with a quark core and neutron 
matter crust in fig. 3b, taking the central energy density as 3200 MeV/fm^. The discontinu- 
ity in the energy density profile for the hybrid star at radius around 3.12 kms corresponds 
to the phase transition from quark matter to neutron matter. As the quark matter EOS is 
softer, it leads to a smaller mass and radius of the quark star. 

For considering the rotational properties of the star, for all the EOS we take the rotating 



star to have the critical value for the angular velocity as given approximately by Qcr = y 
For central energy density of 1600MeV/fm'^, the values of the critical angular velocity are 
1.758,1.551,1.524, 1.254 and 1.667 in units of 10^ sec~^ for Wiringa, Walecka, Bethe- Johnson, 
Sahu et al and Pandharipande EOS respectively. We plot duj/dr as a function of the radial 
distance from the centre of the star with the neutron matter equations of state in fig.4a 
and for quark and hybrid stars in fig.4b. It starts with zero at the centre as given by the 
boundary condition and increases to a maximum and then decreases for all EOS except for 
quark matter EOS. The absence of a peak in the case of quark star is due to the fact that 
there is no crust for such stars and the energy density is finite at the surface even if the 
pressure vanishes, because of the nonzero nature of the bag constant. The height of the 
peak and its location depend both on the crustal thickness as well as the value of the critical 



angular velocity, Vtcr- As has been shown by Datta et al ||2^, the crustal thickness is larger 
for softer equation of state. The height of the peak increases with the stiffness of the EOS 
and increase in the value of VL^r- Wiringa EOS which is the stiffest (for the ec considered) and 
has the highest value for VLcr-, exhibits the maximum peak. For the EOS given by Sahu et 
al, which though is not the softest EOS among those we consider, the star has the minimum 
value for VLcr-, which reflects in the peak being minimum. f2cr for the Walecka model is less 
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than that for Pandharipande EOS, but is stiffer than the latter. These two effects compete 
with each other to make the peaks almost of the same height, where as the position of the 
peak from the surface increases with the crustal thickness as well as with flcr- 

A. Magnetic field topology 

We then consider the field lines for different configurations and study their topology. In 
Fig. 5, for the central energy density, ec=1600 MeV/fm^, we have plotted the magnetic 
field fines for different EOS- (a) Bethe- Johnson, (b) Sahu et al, (c) Walecka and (d) Wiringa 
EOS. The corresponding values for the gravitational potential (GMq/R) are 0.27, 0.284, 0.31 
and 0.33 respectively. The difference in the field line densities seem to indicate that with 
the equation of state getting stiffer, the field line density increases indicating the necessity 
of higher fields to bind the matter configuration. This is consistent with the magnetohydro- 
dynamics assumption of taking into account the concept of frozen in field lines. 

We plot in fig. 6 the magnetic field strength (for ^ = 7r/4) as a function of the radial 
distance in the interior of the star. We have chosen the value for fi approximately as 
1.73x10^^ kms, corresponding to an equatorial magnetic field strength of around 10^^ Gauss 
for a star of radius 10 kms. As can be seen from the graph, the field strength at around 
a distance of 0.1 kms from the center is increased by about six orders of magnitude of the 
same at the surface for all the EOS. 

Rotation also modifies the field strength. We next plot the percentage change in the 
field strength (for 9 — n/A) as compared to non rotating configurations for different EOS. 
As is clear from the Fig 7a, this increase as compared to the nonrotating case is more for 
the stiffer EOS. The field strength can get increased up to 33% due to rotation as one 
approaches towards the core for Wiringa EOS, where as near the surface the increment is 
up to about 10% only. 

We thus have considered here the effect of different EOS on the gross structural properties 
of the rotating neutron stars. We have also considered the case of hybrid stars with a quark 
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matter core and neutron matter crust which have larger mass as compared to the pure 
quark star. Because of the quark hadron phase transition inside such hybrid stars, there is 
a discontinuity in the energy density profile e(r) as well as in the rotational perturbations 
mo(r), po{r), u{r) and duj/dr. However, the discontinuity is more prominent in the energy 
density profile as compared to the rotational perturbation functions. We have also calculated 
the dipole magnetic field structure inside such compact objects. Stiffer EOS yield more 
modification in the magnetic field strength as compared to the nonrotating stars. This is a 
refiection of the fact that the rotational perturbations are more for the stiffer EOS. It may 
be noted that despite the fact that the EOS for Walecka model as well as for Sahu et al 
are in the unstable region for the central energy density chosen for the present calculations, 
the stellar properties as considered here show a consistent pattern. However, it is seen that 
the quark hadron phase transition does not show any discontinuity in the magnetic field 
strength for the hybrid stars. 

To summarize, we have studied rotating compact stars endowed with magnetic fields and 
have discussed some of their gross properties for various EOS including the possibility of a 
hybrid star. However, in the present work, we have not included the effect of magnetic fields 
on the EOS. Such effects with a constant magnetic fields have been shown to soften the EOS 



25| leading to a less massive star. For a more realistic dipole structure of the magnetic field. 



one does not yet know the EOS. 
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Figures Captions 

Fig. 1. Plot of energy density (e) vs. pressure (p) each in units of MeV/fm^. 

Fig. 2. The pressure in MeV/fm^ is plotted as function of baryon chemical potential, hb- 
The sohd line corresponds to the bag model EOS and the dashed line to the Bethe Johnson 
equation of state. 

Fig. 3a. The energy density in MeV/fm^ is plotted as function of the radial distance from 
the centre in kms with the neutron matter EOS. 

Fig. 3b. The energy density in MeV/fm^ is plotted as function of the radial distance from 
the centre in kms for the quark star and the hybrid star. 

Fig. 4a. ^ in 10^/sec/kms vs the radial distance, r in kms with the neutron matter EOS. 
dr 

Fig. 4b. ^ in 10^/sec/kms vs the radial distance, r in kms for the quark star and hybrid 
star. 

Fig. 5. Field lines are plotted in X-Z plane for EOS given by Bethe- Johnson, Sahu et al, 
Walecka and Wiringa in (a), (b), (c) and (d) respectively. 

Fig. 6. The logarithm of the magnetic field strength, B [6 — n/A) in Gauss in the interior 
of the star. 

Fig. 7a. The percentage increase in the field strength as a function of the radial distance in 
kilometers due to rotation of the star with neutron matter EOS. 

Fig. 7b. The percentage increase in the field strength as a function of the radial distance in 
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kilometers for the quark star and the hybrid star. 
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Fig. 6 
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